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ABSTRACT 

We develop a general theory of buoyancy instabilities in the electron-ion plasma with 
the electron heat flux based not upon MHD equations, but using a multicomponent plasma 
approach in which the momentum equation is solved for each species. We investigate 
the geometry in which the background magnetic field is perpendicular to the gravity 
and stratification. General expressions for the perturbed velocities are given without 
any simplifications. Collisions between electrons and ions are taken into account in the 
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momentum equations in a general form, permitting us to consider both weakly and strongly 
collisional objects. However, the electron heat flux is assumed to be directed along the 
magnetic field that implies a weakly collisional case. Using simplifications justified for an 
investigation of buoyancy instabilities with the electron thermal flux, we derive simple 
dispersion relations both for collisionless and collisional cases for arbitrary directions of the 
wave vector. The collisionless dispersion relation considerably differs from that obtained 
in the MHD framework and is similar to the Schwarzschild's criterion. This difference is 
connected with simplified assumptions used in the MHD analysis of buoyancy instabilities 
and with the role of the longitudinal electric field perturbation which is not captured by the 
ideal MHD equations. The results obtained can be applied to clusters of galaxies and other 
astrophysical objects. 

Subject headings: convection - instabilities - magnetic fields - plasmas - waves 
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1. INTRODUCTION 

Physical processes taking place in astrophysical objects are defined by the physical 
parameters of the latter. In many cases, evolution of these parameters can lead to 
instabilities influencing the dynamics of these objects. Convective or buoyancy instabilities 
arising as a result of stratification of astrophysical objects are among those instabilities 
that may operate under different conditions from stellar interiors (e.g., Schwarzschild 1958), 
accretion disks (Balbus 2000, 2001), and neutron stars (Chang & Quataert 2009) to hot 
accretion flows (e.g., Narayan et al. 2000, 2002) and even galaxy clusters and intracluster 
medium (ICM) (e.g., Quataert 2008; Sharma et al. 2009; Ren et al. 2009). Analogous 
instabilities also exist in the neutral atmosphere of the Earth and ocean (Gossard & Hooke 
1975; Pedlosky 1982). It is believed that convective instabilities have a vital role not only 
in transporting energy but in driving turbulence in many astrophysical systems. When the 
entropy increases in the direction of gravity, a thermally stratified fluid becomes buoyantly 
unstable according to the Schwarzschild criterion (Schwarzschild 1958). This well-known 
condition is successfully applied in the theoretical modeling of stellar structure and is valid 
irrespective of existence of the magnetic field. But a stellar fluid is a strong collisional 
system. At the same time, there are systems in which plasma is tenuous and hot such as, 
for example, ICM (Sarazin 1988). Majority of the mass of a galaxy cluster is in the dark 
matter. However, around 1/6 of its mass consists of hot, magnetized, and low density 
plasma known as ICM. The ICM is classified as a weakly collisional plasma with the 
electron number density n e ~ 10~ 2 — 10" 1 cm -3 , the electron temperature T e of the order 
of a few keV (e.g., Fabian et al. 2006; Sanders et al. 2010), and the magnetic field strength 
B ~ 0.1 - 10 //G (Carilli & Taylor 2002). Thus in ICM, the mean free path of ions and 
electrons is much larger than their Larmor radius. Cosmic rays play also an important 
role in the physics of ISM. Recent studies show that centrally concentrated cosmic rays 
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have a destabilizing effect on the convection in ICM (Chandran & Dennis 2006; Rasera and 
Chandran 2008). 

In the recent past, one has included the anisotropic heat flux in weakly collisional 
plasmas and obtained additional instabilities. These instabilities have been shown to arise 
when the temperature increases in the direction of gravity, if the background thermal 
flux is absent (the magnetothermal instability (MTI)) (Balbus 2000, 2001), and when the 
temperature decreases along the gravity at the presence of the latter (the heat buoyancy 
instability (HBI)) (Quataert 2008). The growth rates of MTI and HBI are of the same 
order of magnitude as the growth rates without heat flux. 

Previous theoretical models applied for study of buoyancy instabilities are based on 
the ideal magnetohydrodynamic (MHD) equations (Balbus 2000, 2001; Quataert 2008; 
Chang & Quataert 2009; Ren et al. 2009). However, the ideal MHD does not capture some 
important effects. One such effect is the nonzero longitudinal electric field perturbation 
along the magnetic field. We show here that the contribution of currents due to this small 
(at the present case) field to the dispersion relation can be of the same order of magnitude 
as that due to the electric field components transverse to the magnetic field. Besides, the 
MHD equations can not take into account the existence of various charged and neutral 
species with the different masses and electric charges and collisions of different species with 
each other. On the contrary, the plasma E-approach deals with dynamical equations for 
each species. Starting with Faraday's and Ampere's laws, one obtains equations for the 
electric field components. Such an approach allows us to follow the movement and change 
of parameters of each species separately and to obtain rigorous conditions of consideration 
and corresponding physical consequences in specific cases. This approach permits us to 
include various species of ions and dust grains having different charges and masses. In this 
way, streaming instabilities of rotating multicomponent objects with different background 
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velocities of species (accretion disks, molecular clouds and so on) have been investigated 
by Nekrasov (2008, 2009 a, b, c, d). It has been shown that these instabilities have the 
growth rates much larger than that of the magnetorotational instability (Balbus 1991). 
In some cases, the standard methods used in MHD leads to conclusions that significantly 
differ from those obtained by the method using the electric field perturbations. One of 
such examples concerning the contribution of collisions in the electron-ion plasma has been 
shown in (Nekrasov 2009 c). 

A study of buoyancy instabilities with the electron heat flux by the multicomponent 
E-approach has been performed by Nekrasov and Shadmehri (2010). A geometry has 
been considered in which the gravity, stratification, and background magnetic field are all 
directed along the one axis. Solution of the dispersion relation has shown that the thermal 
flux has a stabilizing effect. The same problem solved by the ideal MHD leads to another 
result (Quataert 2008). 

In this paper, we apply the multicomponent E-approach to study buoyancy instabilities 
in magnetized nonuniform electron-ion astrophysical plasmas in which the background 
magnetic field is perpendicular to the gravitational field. The inhomogeneity is assumed 
to be directed along the gravity. We include collisions between electrons and ions and 
electron thermal conductivity. The dynamical frequency and the wave phase velocity 
along the magnetic field are assumed much less than the ion cyclotron frequency and 
electron thermal velocity, respectively. In this case, the electrons can be regarded as 
inertialess. The pressure anisotropy of species is not taken into account. Thus we exclude 
the pressure-anisotropy-driven firehose (e.g., Vedenov & Sagdeev 1958; Kennel & Sagdeev 
1967 a, b; Schekochihin et al. 2008) and mirror (e.g., Vedenov & Sagdeev 1959, Southwood 
& Kivelson 1993; Califano et al. 2008) mode instabilities. We also neglect the viscous and 
finite Larmor radius effects. The corresponding conditions can be easily obtained from 
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momentum equations given, for example, in (Nekrasov 2009 c, d). Solutions for perturbed 
velocities of species are found in a general form. Using conditions suitable for buoyancy 
perturbations, we obtain simplified expressions for perturbed values. The final dispersion 
relation is derived for the case in which the anisotropic thermal flux is important. The 
growth rates are found for collisionless as well as collisional cases. Solutions of the dispersion 
relation are discussed. 

The paper is organized as follows. In Section 2, the fundamental equations are given. 
The equilibrium state is considered in Section 3. The perturbed ion velocity, number 
density, and thermal pressure are obtained in Section 4. In Section 5, we consider the 
perturbed velocity and temperature for electrons. Perturbed current components are 
calculated in Section 6. In Section 7, we consider conditions to simplify the contribution 
from collisions. The dispersion relation is derived in Section 8. In Section 9, the main 
points of our analysis are discussed. Implications of results obtained for galaxy clusters are 
considered in Section 10. Concluding remarks are given in Section 11. 



2. BASIC EQUATIONS 



We start with the following equations for ions: 



a r + v,-vv,= 



m^rii 



rrii 



rrijC 





the momentum equation, 




(2) 



the continuity equation, and 



— } + Vi • Vpi + 7PiV • Vi = 0, 



(3) 



the pressure equation. The corresponding equations for the inertialess electrons are: 



= -— +g e E+\xB-m e i/ ei (v { -v ! ), (4) 

n e c 

_£ + V ■ n e v e = 0, (5) 

^ + v e • Vp e + 7 PeV • v e = - ( 7 - 1) V • q e , (6) 

&T 1 

+ v e • VT e + ( 7 - 1) T e V . Ve = - (7 - 1) - V • q e , (7) 

at n e 

where the last equation is the temperature equation and q e is the electron heat flux 
(Braginskii 1965). In Equations (l)-(7), qj and rrij are the charge and mass of species 
j — i, e, Vj is the hydrodynamic velocity, rij is the number density, pj = nfTj is the thermal 
pressure, Tj is the temperature, v ie (v ei ) is the collision frequency of ions (electrons) 
with electrons (ions), g is the gravity, E and B are the electric and magnetic fields, c is 
the speed of light in vacuum, and 7 is the adiabatic constant. As for the electron heat 
flux, we consider the case of weakly collisional plasma for which the electron cyclotron 
frequency u ce = q e B/m e c is much larger than the electron-electron collision frequency v ee , 
i.e. u ce ^> is ee . In this case, the electron thermal flux is mainly directed along the magnetic 
field, 

q e = - Xe b(b-V)T e , (8) 

where Xe is the electron thermal conductivity coefficient and b = B/_B is the unit vector 
along the magnetic field. We only take into account the electron thermal flux (8) because 
the ion thermal conductivity is considerably smaller (Braginskii 1965). We also assume that 
the thermal flux in equilibrium is negligible. 

Electromagnetic equations are Faraday's 
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and Ampere's 




(10) 



laws, where j = ^ qjUjVj. We consider wave processes with typical time-scales much larger 
than the time the light spends to cover the wavelength of perturbations. In this case, 
one can neglect the displacement current in Equation (10) that results in quasi-neutrality 
for both the electromagnetic and purely electrostatic perturbations. The magnetic field 
B includes the background magnetic field B , the magnetic field B 0cur of the background 
electric current (when it presents), and the perturbed magnetic field. 

Our basic equations are the same as that in (Nekrasov & Shadmehri 2010). However, 
in the present case, the equilibrium state differs from the equilibrium in (Nekrasov & 
Shadmehri 2010), where magnetic field lines and the gravity are parallel to each other. 
In the next section, we analyze the equilibrium state in which the magnetic field is 
perpendicular to the direction of gravity. 



At first, we consider the equilibrium state. In this paper, we study the configuration 
in which the background magnetic field B is parallel to the z-axis, and the gravity and 
stratification are parallel to the x-axis. Let, for definiteness, g be g = — xg>, where g > and 
x is the unit vector along the x-direction. Then, Equations (1) and (4) are the following: 



3. EQUILIBRIUM 



3.1. Background velocities 



Vj • Vv i0 




v i0 x B - i/? (v i0 - v e0 ) , 



(11) 



miUio 







Vp e0 



+ q e E + — v e0 x B - m e z/° (v e0 - v i0 ) , 



(12) 
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where the index denotes the background values. Our system is uniform along the y- 
and z-axes. Therefore, we assume that E 0y = E 0z = 0. However, for generality, we keep 
Eq x 7^ 0. We also assume that v i0z = v e0z = 0. Then, from the y-components of Equations 
(11) and (12), we obtain 

Vi0x=Ve0x= (ViOy ~ V e 0y) (13) 

provided that u C i = qiB /miC ^> dv^y/dx. In this case, the contribution of the ion inertia 
along the y-axis produced by the shear of velocity v^ y along the x-axis can be neglected 
in comparison with the y-component of the Lorentz force. From the x-components of 
Equations (11) and (12), we also find under condition u ci Vio y ^> dvf 0x /2dx 

ViOy =v E + v g + v di = -c— H , (14) 

-DO U c i UJ c i 



Eqx Qe 

VeOy = V E + V de = - C— H , (15) 



where 



1 dpio 1 dp e0 
9i = 7T- ,9e = 7T- ■ (16) 

miUio ox rriin e o ox 

The velocities ve, v g , and Vdi, e are the electric, gravitational, and diamagnetic drifts, 
respectively. We assume that = — q e . If values gi and g e are inhomogeneous and gi ~ g e , 
then conditions given above to justify solutions (13) and (14) can be written in the following 
form: 

+ (17) 
where pi = VTi/oj ci {vt% is the ion thermal velocity) is the ion Larmor radius and 
L^ 1 = dg^ e / g^ e dx (an uncertain contribution of the electric drift is not taken into account 
in condition [17]). We note that there are no restrictions on the relation between u C i and 
z/° e in inequality (17), i.e. the case uf e > u ci can also take place. 
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3.2. Continuity and pressure equations 

Having the background velocities, we can now consider the ion continuity equation (2). 
In equilibrium, we have 

~df dx nioVi0x = ^ ' 

This equation is nonstationary because of the x-dependence of v i0x . The typical time 
to of the density evolution is to ~ ^/^ie^Ph where k = drio/nodx (n^o = n e o = no). 
When studying the linear perturbations, we will not take into account that the medium is 
nonstationary. This means that we will consider perturbations which develop much faster 
than t - Thus, we can safely assume that the medium is stationary. The same also relates 
to the ion pressure equation and corresponding equations for electrons. 



4. LINEAR ION PERTURBATIONS 



4.1. Equation for the perturbed ion velocity 

Let us write the equation of motion (1) for ions in the linear approximation, 

dvn Vpu Vp i0 rin q t qi 

— — + v i0 • Vvji + Vii • Vv i0 = + + — Ei + v a x B (19) 

dt miUio miUiQUio rrii rriiC 

Qi 

H ^V i0 xBi- I/? (Va - V e i) - v\ e (V i0 - V e0 ) , 

rriiC 

where the index 1 denotes the values of the first order of magnitude. The diamagnetic drift 
in the stationary velocity Vj je0 is not a real velocity (e.g. Nishikawa & Wakatani 1990). 
However, this drift must be taken into account in the hydrodynamical equations evoking 
drift waves (e.g. Vranjes et al. 2003 for the same initial state). The last collisional term in 
Equation (19) appears as a result of perturbation of collision frequency due to the density 
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and temperature perturbations: v\Jv® e = n e \/n e o — 3T e i/2T e0 . This effect due to the density 
perturbation has been involved by Nekrasov (e.g. 2009 b). 

Below, we do not include the shear of the background velocity on the left hand-side 
of Equation (19) (the same also relates to the continuity and pressure equations). It is 
possible, if 

A* = ^ + v l0 -V»— . (20) 

The first inequality (20), D ti ^> dvio x /dx, coincides with the condition of medium 
stationarity (see Eq. [18]). The second inequality, D ti ^> dv^y/dx, is satisfied automatically 
ifvioy = const. When v^ y ^ const, this condition can be written in the form D ti ^> uj C ipf/L 2 . 
We further introduce the following notations: 

F a = A El -^( Vil -v el ), (21) 
rrii 

Gn = Fji + — ^-Vj x Bi, 

rriiC 

Gil = v} e (Vj - v e0 ) . 

Then Equation (19) takes the form 

AiV 41 = -^L + ^^+G tl - C* + -^-v a x B . (22) 
miUio miUio n i0 rriiC 



4.2. Perturbed ion continuity and pressure equations 

In the linear approximation, the ion continuity (Eq. [2]) and pressure (Eq. 
equations are given by 



Dunn + vnJ^- + n i0 V • v a = 0, (23) 
ox 
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DtiPn + v ax-^r + 7PioV • v a = 0. (24) 

Equations (22)- (24) are used to find perturbed ion velocity components. These calculations 
are given in the Appendix A. For simplicity, we adopt that g i>e = const (see Eq. [16]). 
This is true when the gravity is directed along the background magnetic field (Nekrasov & 
Shadmehri 2010). This case is also satisfied, if we impose the condition v^ ely = (see Eqs. 
[14] and [15]). However, in a general case g iye ^ const. 

5. LINEAR ELECTRON PERTURBATIONS 
5.1. Equation for the perturbed electron velocity 

Consider Equation (4) in the linear approximation 



Vp e i Vp e0 n el q e q e l 

= 1 hg e Ei + — v e i xB H — v e0 xB 1 -m e v ei (v el - vn)-m e v ei (v e0 - v i0 ) . 

n e o n e0 n e0 c c 

(25) 

Here z^/V^ = nn/riio — 3T el /2T e0 . We introduce the following notations: 

F e i = g e Ei - m e u° ei (v e i - Va), (26) 

G e i = F el +— v e0 x Bx, 

c 

C e i = (V e0 - Vfl)) ■ 



Then Equation (25) takes the form 



= _v^ + v ?M!?fl + Gei _ ai + ^ ixBo 

n e0 n e0 n e0 c 



5.2. Perturbed electron continuity, pressure, and temperature equations 

The continuity equation (5) in the linear approximation is given by 

D te n el + v eixj^- + n e0 V ■ v el = 0, (28) 

where D te — d/dt + v e0 • V. Here, we also neglect the contribution of equilibrium velocity 
inhomogeneity, i.e. we assume that D te ^> dv e0x /dx. 

The linear electron pressure equation (6) is the following: 

DtePel + V ^~^- + 7PeoV • V el + (7 - 1) V • q el = 0, (29) 

where q e i is the linear electron thermal flux. 

The perturbed electron temperature equation (7) has the form 



dT 1 
D te T el + v elx —^ + (7 - 1) T e0 V • v el = - (7 - 1) —V • q el . (30) 
ox n e0 



5.3. Perturbed electron thermal flux 



Expression for the thermal flux (8) in the linear approximation is 



where 



dT 

q e i = -Xeob (b r V) T e0 - Xeob -^-, (31) 



, B 1 B lz 
bi = — b - 



Bo B 
Thus from Equation (31), we have 



q e ix, y = 0, (32) 
dT e Q B\ x dT e \ 
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Using Equation (32), we find the value V • q el as follows 

_ dT e0 1 dB lx d 2 T el 

V-q el = -x e o^^^-Xeo^. (33) 

The first term on the right hand-side of Equation (33) can be considered as a source of 
the heat, and the second term describes the thermal diffusion. This expression for V • q el 
is analogous to the one arising in the case in which the stratification is directed along the 
background magnetic field and there exists an anisotropic background thermal flux (see 
Nekrasov & Shadmehri 2010, Eq. [49]). 



5.4. Perturbed electron temperature and pressure equations with the thermal 

flux 



Equation (30), taking into account expression (33), can be written in the form 

-l 

B\ x , 



(Dte + ft) T e i = -Velx^TT^ ~ (l ~ 1 ) Te ° V ' Vel ~ -^^^ ( ^~ 



dx 



dx B \dz 



where the operator 



1 d 2 



is introduced. 



(34) 



(35) 



Equation (34) permits us to define the value V • q el . Substituting this value in 
Equation (29), we obtain the desirable equation for the perturbed electron pressure 







DtePel + 



dp, 



eO 



ft 



dx 



- n e0 



eO 



(Ae + ft) dx 



J Dte + ft / x 

^ + Pe0 (A e Tft) v ' Vel (36) 



+n e o- 



dT e0 1 A,ft 



dx B (Ae + ft) \dz 



d 



B 



lx ■ 
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Equations (27), (28), and (36) are used in the Appendix B for determining the velocity 
v el . For convenience and compactness of all the equations and expressions, we formally 
apply the method of division by time and spatial derivatives (see Eqs. [34], [36], and below). 
This permits us not to use the Fourier transformation at the initial stage of calculations, but 
rather gives a possibility to obtain the spatiotemporal differential equations for variables 
under study (for example, for numerical calculations). When the dispersion relation is 
derived in Section 8, all derivatives will be changed by their Fourier-images. 



6. CURRENT ji 



Expressions for Vji and v el obtained in the Appendices A and B have a sufficiently 
general form. In particular, one can take into account the background velocities of species 
that results in electromagnetic streaming instabilities. Some studies of such instabilities 
in astrophysical objects are given in Nekrasov (2007; 2008; 2009 a,b,c,d). However, in the 
present paper, we do not consider this effect, assuming that 

D t > v iiC0 -V, (37) 
AEi > Vv i)e0 -Ei, 

where D t = d/dt. Thus, D ti:S ~ D t . The second inequality (37) is obtained from condition 
cEi 3> Vi : eo x Bi (see Eqs. [21] and [26]) and by using Equation (9). Under conditions 
(37), one can also neglect the contribution of Cj je i in Equations (22) and (27). We further 
consider perturbations with the dynamical frequency much smaller than the ion sound 
frequency. Keeping the spatiotemporal derivatives, this means that c 2 si d 2 /dz 2 3> Df . We 
also assume a local approximation when the perturbation wavelength in the direction of 
inhomogeneity is smaller than the inhomogeneity scale, i.e. d/dx ^> g ije /(? si e . Below, we 
find components of the perturbed current ji. For simplicity, we omit the index by v® e . 



- 16 - 



6.1. Current j lx 



From Equation (A15) and under conditions given above, we obtain 



A / d 2 d 2 

Vilx = ~2 I — + 



u li \dy 2 dz 2 J \dz 
u ci dy \dz) llz bj 2 ci 



D u d 2 (d 



u 2 ci dxdy \ dz ' 1 



(38) 



ujrrDt d ( d 



4 d v \ dz 



+ 



d_ 

dx 



d_ 

dz 



ilzt 



where the additional condition d/dx 3> {gi/c 2 si ) (d/dy) 2 (d/dz) 2 has been assumed. From 
Equation (B4), we have 

1 „ 



Velx z 



ely 



1 9( l\ F*, 



(39) 



m e u) ce " m e u ce dy \dz 
Expression (38) coincides with expression (39) in the limit rrii — >■ at substitution i -H- e 
and rriiFn ->• F el . 

Using expressions (38) and (39), we can find the current ji x . It is convenient to find 
the value 4irji x /D t . As a result, we obtain 



4.7T 



D f u, 



■ --, i ~1 I 9 O'xx I jix O'XxEix Qj X yE\y d xz E\ z ~\~ o O'xyjly 9 Q*xzjlz- 



pi 



4:71 Ui, 



(40) 



pi 



The following notations are introduced here: 



wi- / a 2 <9 2 



w 2 j ^<9?/ 2 (9^ 2 J \dz) ,Qxy oj 2 ci dxdy \dz 



d \ 



-2 



coi d 2 (d 



pi 



-2 



(41) 



uJciDt d fd_\ + d_ 
c 2 si dy \dz J dx 



d_ 

dz 



When deriving Equation (40), we have used expressions (21) and (26) for F a and F el , 
respectively. We note that in Equation (40) and below in this section, we have not imposed 
any restrictions on the collision frequency V{ e in comparison with the ion cyclotron frequency 
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6.2. Current j ly 



Let us consider the velocity vn y . Using corresponding conditions, we can write from 
Equation (A18) 



Vily = 



L Fa Ei 92 ( d 

uj 2 ci dxdy \ dz 



ilx 



(42) 



Dt 



d 2 



+ 



d 2 



dx 2 dz 2 1 \dz 



d 



ily 



1 d ( d \~' 1 
+——( — ) F ilz + — 



u ci dx V dz 



Dj d ( d 



c 2 , dx V dz 



9% 



Here, we have assumed that d/dx 3> (gi/c 2 si ) (d/dx) 2 (d/dz)~~ 2 . We keep the last term in 
Equation (42) because the term proportional to (00^/00^) F ilz can be canceled with the 
corresponding electron term (see below). Here and below, we have introduced the notation 

del . 



bi,e 2 

si.e 



9i,e (7 - 1) + 



si,e 



dx 



(43) 



which may be named the ion (electron) Brunt- Vaisala frequency. 



The velocity v ely is defined by Equation (B6). For the case under consideration, we 
obtain 

-1 



7 



1 



elz 



(44) 



2 7 
-rrnu be 



d ( d 



1 

^se 

+g e mi 



^D t + ft m e oj ce dy \ dz 

(7 - 1) 9e 



dc 2 

uu se 



7 A + ^ dx 
dT e0 1 n ( d 



Fell 



d 



dz 



elz 



B 



lx- 



T e0 dx B ~fD t + Q, \dz. 
If we put Q = in Equation (44) and take m ; — > in Equation (42), we obtain the full 
conformity of both equations with each other. 
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Using expressions (42) and (44), we find the current j ly . Taking into account (21) and 
(26), we obtain 



in 



1 + 



DM* 
uj 2 - 



' !/!/ 



Jiy 



—a yx E lx + [a yy + a* yy ) E ly 



(a V z ~ a* yz ) E 



Iz 



(45) 



4-KV,; f 



UJ 



pi 



UJ 



2 a yzjlz; 



pi 



where the following notations are introduced: 



<4 d 2 



l yx 



l yz 



d 



uj 2 ci dxdy \ dz 



i a yy 



uj 

UJ 



d 2 

, dx 2 dz 2 



d 



<i 



<4e 



dz 



UJ 



pi 



D 2 d 



D t u ci I c 2 si dx \dz 



d 



-2 



UJ 



9i 



(7-1) ft 



D 2 D 2 7 D t + n 



dc 2 



(,46) 



se 

d ( d 



7 A + Q dx 




ujIX ^ A 2 D t (jD t + n)S dy\dz 



1 di 



"yy 



p% 

D 2 ujJ t 9e T e0 dx 1 D t + n , ^ yz ~ ™ yy dy \dz. 

The terms proportional to a* y z are connected with the contribution of the heat flux ~ B lx 
(see Eq. [32]). This magnetic field perturbation has been expressed via electric field 
perturbation through Equation (9). 

We see from the expression for a yz that for Q = and T i0 = T e0 the last two terms in 
the first figured brackets are canceled. 



dT 



eO 



d ( d 



6.3. Current j lz 

From Equation (A20), we find the simplified ion velocity vn z 
dv ilz 1 dF ilx D 2 d ( d\ 2 D t dF ilx 1 dF ily g t 



dz uj ci dy uj ci c 2 si dy \dz J uj 2 ci dx uj ci dx 

D 2 d ( d\\ 1 ,_ 2 , *.dF ay D t fd^- 



0J C iC 2 si dx V dz 



( 9 \ f 1 <n? + ^) dFa > Dt ( 9 \ f 
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We have taken into account that d/dx 3> (gi/c 2 si ) [(d/dx) 2 + (d/dy) 2 ] (d/dz) 2 . From 
Equation (B8), we have 



dVeU 

dz 



1 dF elx 1 dR ( 



ely 



+ 



m e uj ce dy WeWce dx C 2 el 



9e + 



dP, 



ely 



D, 



(jD t + £1) \m e u ce J dy c 2 sel rrii \dz 
dT e0 1 D t n 



n dT e0 

(D t + Q) rriidx 
i 

Relz 



m e u ce 



ely 



(48) 



d 



d 



where 



c 2 sel rriidx B (D t + ft) \dz ' 

2 _ 2 (-yDte + Q) 



(49) 



If we take in Equation (47) rrii — > and put f2 = in Equation (48), we obtain the 
conformity of these two equations. 

Using Equations (9), (47), and (48), we find the longitudinal current j\ z 

4vr / D t v ie \ . . m . 4irv ie . 

77 1 —a zz )ji z = a zx E lx + [a zy + a ) E ly - (a zz + a zz ) E lz —a zx j lx [bO) 



pi 



pi 



2 a zy3ly 



pi 



The following notations are introduced here: 



O'zy 



Si 

UJ 2 - 
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u ci D t d ( d 



4 d v \ dz 



9i_ 

UriDt 1 <?■ 



,2 



c 

si ^sel 



d 




dx 


(dz) 
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dT e0 



(D t + ft) rriidx 



^si ^sel 



D 2 d ( d 



U) 



pi 



" 2 a D 2 



.2 2 2 'jDt \ d ( d 

)t+UJ - + ^WTn)dy-{d-z 



c 2 si dx V dz 
i 



d_ 

dz 



UJ. 



a zy — 



pi 



dT. 



eO 



n 



d 



,a 



d ( d 



u ci D t c 2 el midx (D t + Q) \dzj ' ~ zz a * y dy \dz / 

The terms proportional to a* y in Equation (50) are connected with the contribution of the 
electron current induced by B lx (see Eq. [48]). 
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7. SIMPLIFIED COLLISION CONTRIBUTION 



To take into account collisions between ions and electrons, we involve some simplifying 
assumptions. First of all, we assume that 

1 » -V £ - (52) 

This condition allows us to neglect the collision contribution on the left hand-side of 
Equations (40) and (45). We note that the relation between u ci and z/j e in inequality (52) 
can be arbitrary because D t <C uj ci . For estimations of contribution of the different terms, 
we assume here and below that D t > uj bi , e - We also consider that 

max J (i y 2 . £* l \ L » (i) ~ 2 (53) 

I c 2 si dy \dz) 'ucidxj u 2 ci dxdy \dz J 

where L is the typical inhomogeneity scale length along the x-axis. Under conditions (52) 
and (53), we can neglect the contribution of the collisional terms proportional to uf e . Then 
solving Equations (40), (45), and (50), we can express the current ji through Ei 

47T 

~j^~jlx ^xxEix &xyE\y E xz Ei z , (54) 

47T . 

~T7~3ly &yxE\x &yyE\y + Ey Z ±Li\ z , 

Dt 

47T . 

~p7~jlz E-zxE\x ~\~ ^zyEiy £ ZZ E\ Z . 

The following notations are introduced here: 

Dt Vie «^ K> + <y) 

&XX ^xxi&xy u-xy 2 /1 -J \ ' V / 

ui pi {1 — ia z ) 

_ D t v ie 

&XZ &XZ 



1 1 2 



n (n - n* \ — i^zz + Q* 2 ) 

a xy {a yz a yz ) (1 _ i( y 



D t u ie a yz a zx _ , * _ D t Vie a yz {a>zy + a* zy ) 

£yx ayx+ w 2 (i- idz y £ yy- a yy + a yy w 2 (i_ i( y 

„ gfe (a zz + a* J 
£yz - a yz a yz + ^2 (l-id,) ' 



id z ' £zy 1 — 2<i z ' £zz 1 — i<i 2 
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where 

D u° 

id z = a zz -^. (56) 



We keep in = x,y,z) some terms, which will turn out to be small, because we do 

not know a priori solution of the dispersion relation. 



8. DISPERSION RELATION 



To derive the dispersion relation, we apply the Fourier transformation to the 
electromagnetic field and current, assuming that Ei~ Eifc exp (ikr— iout), where k = {k, u}. 
We take into account all three components of the wave vector, k = (k X) k y , k z ). Then using 
Equations (9), (10), and (54), we obtain 



AfcEi fc — 0, 



where E lfc = (Ei xk , E lyk , E lzk ) and the matrix A k has the form 



(57) 



A fc = 



n 2 y + n 2 z - e a 



~Tl x Tly ~\~ &xyi ^x^z ~t~ ^xz 



-yyi 



^x^z &zxi HyTl z &zyi 



r n>y r n> z ^yz 

£ zz 



n 2 x + n 2 y 



Here n = kc/cu. We keep for values the same notations as above. The dispersion relation 
is the determinant of the matrix A k equal to zero 



(jly 4~ fl z £ xx j [(^a; 4" H z &yy) &zz £zy&yz\ 

(jlxTly ^yx) [(jlxTly ^xy) ^zz ^xz^zy ^zx^-yz] 



(58) 



When obtaining Equation (58), we have taken into account that according to Equations 
(41), (46), (51), and (55) inequalities e zz > e xx ,e xy ,e xz ,e yy ,e yz ,e Z y and e xx e zz > e xz e zx are 
satisfied. Expression e xz e zy — e zx e yz can be omitted under the condition 1 > 1/(1 - id z ) kL, 



-22- 



when > w or Tj ^ T e at 11 C w. Here, the value k is equal to k — |k| and 
n = ( 1 -l) Xe0 k 2 z /n e o (seeEq. [35]). 

To further simplify the dispersion relation (58), we neglect the collision contribution in 
e xy and e yx . The corresponding condition can be written in the form 



1 > 



Vie \ I 1 W 



kL k"lc 2 si 



(59) 



(1 - id z ) 

We remind that the ratio u ie /oj C i in this inequality is arbitrary. Then Equation (58) takes 



the form 



ft z C A + , ,2 "^t I £ 3/3/l ^ 



£ zy^-yz 



y (A£ + *jj + *2)' 



(60) 



1/2 y 

where = (Bq/ Aitrriinio) is the ion Alfven velocity. The following notation is introduced 
here: 



'3/3/1 



+ 



be 



tA 



A 2 -D 2 7 A + 



°3/3/- 



A^ie <V ( a zy + a *zy) 



Ul. 



pi 



id. 



(61) 



For convenience, we retain the symbol D t = —iu on the right hand-side of Equation (60). 
Below, we consider the dispersion relation in the collisionless as well as in the collisional 

cases. 



8.1. Collisionless case 

In this case, we assume that 

d z < 1, (62) 

where d z is defined by expression (56). Under condition (62), the collisional term in e yyl is 
unimportant. Then the dispersion relation is given by 
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where c 2 = c 2 si c 2 sel / (c 2 si + c 2 sel ) and 



W = ^ I w 2 + ul 



7A 



+ <7, 



<9T 



eO 



ft 



6i 



< -— ^ 

ft C se 



7 A + ft T e0 dx-fD t + Q 

dc 2 



The condition 



has been supposed. 



l > 



Wci Lk 2 (D t + n) 



9i_ 

r2 



ft 



(64) 



The case ft S> A In this limit, the value W is the following: 



u bi + ft 



dT 



eO 



T e0 &r 



ft 



1 dc 2 



c 2 se dx 



This expression can be rewritten in the form 

(ft + ft) 



W = 7- 



r 2 _2 



( 7 -l)^+ 7 



ft_ 

^2 



7ft 



eO 



(65) 



where we have used expression (43) for cu 2 i e . We see that instability, W < 0, is possible, if 
the temperature increases along the gravity (we assume that g i;S has the same sign as g). 



The case ft <C D t In the case O C Dj, the thermal conductivity is absent. The value W 
is given by 



or 



(ft + ft) 

r 2 r 2 



7 /<9T j0 , &r t 



(7-l)(ft + ft) + - + 



eO 



(66) 



Comparing Equations (65) and (66), we see that the thermal conductivity is not of 
fundamental importance. 
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8.2. Collisional case 

We now assume that 

id z > 1. (67) 

In this case, we can neglect the term e zy e yz / e zz in Equation (60). However, the collisional 
term in expression (61) gives the same contribution as other terms. As a result, we obtain 
again Equation (63) with W defined by Equation (64). Thus, the dispersion relation is 
the same for both the collisionless and collisional cases. We note that this result has also 
been obtained for the case in which gravity is parallel to the magnetic field (Nekrasov & 
Shadmehri 2010). 



8.3. Polarization of the electric field perturbation 



The dispersion relation (58) without e zx e yz —e xz e zy can be obtained, if we neglect some 
terms in the matrix A*. 



A, 



TLy -\- Tl z &xxi ^x^y &xy> 

^rix^iy ~t~ £yxi W'x ^yyi ^yz 



^zyi 



Then Equation (57) takes the form 



(n 2 y + n 2 z - e xx ) E lxk + 
(—n x n y + e yx ) E lxk + (n x + n 2 — e yy ) E lyk — e yz E lzk 

&zyE\ yk ~\- B zz E\ zk 



0, 

o, 

0. 



(68) 



If we put n x = 0, then the magnetosonic wave u 2 = k 2 c 2 A is split, and Equation (60) 
describes the Alfven type wave. When n y = 0, the Alfven wave ui 2 = k 2 c\ is split, 
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and Equation (60) describes the magnetosonic type wave. In the case n x = n y = 0, 
both waves are of the Alfven type. In a general case, we have u 2 ^ k 2 c\ and E lxk = 
[k x k y / (k 2 y + k 2 z )] E lyk . 

We see from the system of equations (68) that E\ zk = ie zy /e zz ) E\ yk <C E\ yk . However, 
the contribution of the longitudinal electric field E lzk must be taken into account in the 
collisionless as well as in the collisional cases (see Eqs. [60] and [61]). 



9. DISCUSSION 



The dispersion relation (60) emphasizes the important role of the perturbed current and 
electric field along the background magnetic field. These perturbations produce the terms 
connected with e yz , e zy , and e zz in the dispersion. Taking into account these terms allows 
us to derive the correct dispersion relation which has the form (63). The value W , defined 
by expression (64), is available for both cases when the electron thermal conductivity is 
present or absent. Expressions (65) and (66) show that the thermal conductivity is not of 
fundamental importance for the buoyancy instability, if the latter can be excited. For the 
instability, the temperature gradient of ions and electrons must have the sign opposite to 
that of <7j i6 . Under assumption that and g e have the same sign as g, the temperature 
must increase along the gravity for the instability to be excited. 

The dispersion relation (63) takes into account collisions. Except for the anisotropic 
thermal conductivity adopted in this paper, the relation between u ci and v ie can be arbitrary 
in the framework of conditions (52), (53), and (59). Under these conditions, the dispersion 
relation for the collisionless (62) and collisional (67) cases is the same. This result has also 
been obtained for the case in which the background magnetic field and gravity are parallel 
to each other (Nekrasov & Shadmehri 2010). 
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From Equations (23) and (24), we can find the ion number density and pressure 
perturbations. Using Equations (38), (42), and (47), we can calculate the value V • v^. 
Keeping the main terms, we have 

V • v a ~ -AfFiiy - ^ F ilz , (69) 

Substituting v ilx and V • Vji in Equation (23), we obtain an estimation 

rin 1/1 driin Qi \ D t ( d \ 1 , . 

- A— ^ — I ^ + ^ J F ily r ( 77" ) 70 

n i0 uJd \n i0 ox ciJ c 2 si \dz J 

We see from Equations (69) and (70) that V- v,i ~ D t nn/niQ. From the system of equations 
(68), it is followed that E lzk = (e zy /e zz ) E lyk . Thus, both terms on the right hand-side of 
equation (70) are of the same order. An estimation for the ion pressure perturbation is 
given by 

J - / t — 9 I o I 1 llz- 

IPiO c S i \ dz J 

We see that na/riio ~ pn/pio due to the longitudinal electric field perturbation. The 
analogous conclusion can also be made for electrons. This result contradicts an assumption 
that rii/rio ^> pi/po which one uses in the MHD analysis of buoyancy instabilities. We note 
that the ideal MHD does not involve the field E z . 

Let us discuss the relevance of conditions used in this paper to real astrophysical 
systems. As an example, we will consider an intracluster medium. However, our conditions 
have a more general applicability. Observations show that all astrophysical objects in 
cosmic space have magnetic fields of /j,G strength in galaxy clusters and molecular clouds 
(e.g., Carilli & Taylor 2002; Wardle & Ng 1999) to G and more in acretion disks (e.g., 
Desch 2004; Donati et al. 2005). For such magnetic fields, the ion Larmor radius p; is 
considerably smaller than the typical inhomogeneity length L in these objects, i.e. L. 
For example, if we take for ICM Bq ~ 1 pG and Tj ~ 1 keV, we obtain p« ~ 2 x 10 4 km. 
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The magnitude of L is ~tens of kpc. In ICM, the collision frequency v ie is much less than 
the ion cyclotron frequency u ci , u ie <C u ci . In dense molecular clouds and accretion disks, 
the ions can be unmagnetized, u ie ^> u ci (Wardle & Ng 1999). However, many conditions 
of our consideration (see below) can also be satisfied in the last case. For example, this 
relates to inequality (17) justifying solutions (13) and (14) in the equilibrium state. The 
electromagnetic buoyancy perturbations have a dynamical frequency u and wavelength A 
much less than the ion cyclotron and sound frequencies and inhomogeneity scale length, 
respectively. Thus, the conditions (52) and (59) are satisfied for a weakly collisional plasmas 
as ICM and can be satisfied for a strong collisional plasma when v ie > u ci . Accounting 
for the fact that \u\ ~ g/c s (indices i and e are omitted), we can treat the medium as 
a stationary one (see Eq. [18]), if 1 3> (is ie /u ci ) (pi/L) that is justified for astrophysical 
objects. The condition (20) is satisfied for the last inequality (which also relates to Equation 
(28)) and under p t <C L in the case v^y ^ const. The condition (53) can be written in the 
form 1 ^> {v ie /u) C i) (X/L; pi/ X), where A < L and pi <C X. The first condition (37) is the 
following: 1 ^> (pi/X) (1; u ie /u C i). It is also true for the transverse component of the second 
condition (37). However, taking into account that E lzk = (e zy /e zz ) E lyk (see Sec. 8.3), 
both sides of the ^-component of the second inequality (37) can be of the same order, if 
v i,eo 7^ 0. Owing to indefiniteness of the background velocities also containing the electric 
field, we here do not consider their effect. Nevertheless, the streaming instabilities can also 
take place. 

In our analysis, we consider for generality that ions and electrons have the different 
temperatures. However, in Eqs. (3), (6), and (7), the terms describing the energy exchange 
between species due to their collisions have not been taken into account. This is possible, if 
the dynamical time scale is smaller than the time scale of smoothing of the ion and electron 
temperatures, i.e. u ie <^ u>. In the opposite case, u ie 3> u>, the perturbed temperatures of 
electrons and ions are almost equal each other. Equations (6) and (7) for electrons will keep 
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their form because v e i Vji. In the case T e0 ~ T i0 , these equations will stay the same with 
the heat flux two times less than the former one. Equation for the ion temperature will not 
be needed. 

10. IMPLICATIONS OF THE OBTAINED RESULTS FOR GALAXY 

CLUSTERS 

According to a simplified point of view, the core of cluster must be cooled. In fact, 
observations show that in most clusters the cooling time-scale near the center of cluster, 
10 8 to 10 9 yr, is much shorter than a cluster's age 10 10 yr (e.g., Fabian 1994; Peres et 
al. 1998; Allen 2000). But such high rates of cooling has been definitely ruled out by 
the X-ray observations showing that the cluster core is sufficiently hot (e.g., Allen 2000). 
Thus, some heating mechanisms are in operation, though we have little knowledge about 
them. Different heating mechanisms from AGN feedback to cosmic rays and turbulence are 
proposed to resolve the "cooling flow problem" (e.g., Eilek 2004; Binney & Tabor 1995; 
Loewenstein et al. 1991; Reynolds 2002). Exactly for this reason, a study of buoyancy 
instabilities in ICM has a purpose to find a solution of this longstanding problem in clusters 
of galaxies. The plasma in ICM is turbulent (e.g., Loewenstein & Fabian 1990; Cattaneo 
& Teyssier 2007). Buoyancy instabilities could be one of possible sources of turbulence 
resulting in the emergence of heat fluxes along the magnetic field in the direction of core. 
Thus, the latter could be heated. The total nonlinear picture of this process which can 
include reorientation of the magnetic field (Parrish & Stone 2007) is very complex for the 
analytical consideration and can only be investigated numerically. 

In the framework of the linear ideal MHD, new buoyancy instabilities have been found, 
when the thermal conduction is the dominant mode of heat transport (Balbus 2000, 2001; 
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Quataert 2008). In the present paper, in the framework of the multicomponent E-approach, 
we have found generalized growth rates (for the same geometry as considered by Balbus 
2000) for both the negligible and dominant thermal conduction. In both cases, the growth 
rates have the same order of magnitude. However, our modified conditions of instability are 
different from the MHD case because of the multifluid nature of the system (see Section 
9). We have shown that conditions for the buoyancy instability have the form which is 
analogous to the Schwarzschild criterion (Schwarzschild 1958). Thus, the ICM plasma can 
be buoyantly unstable for both large and small thermal conductivity. This increases the 
range of wavelengths of the unstable perturbations because the thermal flux is proportional 
to the wave number. Our linear analysis shows that the multifluid nature of plasma in 
the ICM and other astrophysical objects can not be neglected when one investigates heat 
flows as a result of buoyancy instabilities. We also note that buoyancy instabilities are of 
the electromagnetic nature. Therefore, they can contribute to the magnetic field activity in 
astrophysical objects. 



11. CONCLUSION 



In this paper, we have investigated buoyancy instabilities in the magnetized electron-ion 
astrophysical plasmas in which the background magnetic field and gravity are perpendicular 
to each other. We have applied the multicomponent E-approach in which the dynamical 
equations for ions and electrons are solved separately via the electric field perturbations. 
The perturbed current and Faraday's and Ampere's laws have been used to derive the 
dispersion relation. We have included collisions between electrons and ions. Except for the 
anisotropy of the electron heat flux adopted in this paper, in other respects, the relation 
between the ion-electron collision frequency and ion cyclotron frequency can be arbitrary 
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in the framework of approximations, which have been made. The important role of the 
longitudinal electric field perturbations, which are not captured by the MHD equations, 
has been shown. The obtained growth rates for cases of strong and weak electron thermal 
conductivity show that an instability is possible when the temperature gradients of ions and 
electrons are directed along the gravity. We have shown that the relative perturbations of 
number density and pressure are of the same order as a result of action of the longitudinal 
electric field perturbation. 

Results obtained in this paper are applicable to the magnetized collisional stratified 
objects and can be useful for a search of sources of turbulent transport of energy and matter 
in the ICM and other astrophysical objects. 
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APPENDIX A 



A. SOLUTION OF EQUATION (22) 

The components of Equation (22) are the following: 

1 dpn rin ( , . 

L>uv ilx = gi h G ilx + UdVuy, (Al) 

mifiio Ox n i0 

DtiVily— h Gily — Cily ~ W C j^jlX, (A2) 

niiUiQ ay 

D u v ilz = — ^- + G ilz . (A3) 
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Let us apply the operator D ti to Equation (Al). We do not differentiate the value l/n i0 and 
change D ti dpn/dx by dD u pn/dx according to condition (20). Then we use Equations (A2), 
(23) and (24). As a result, we obtain the following equation connecting v ilx and V • v^: 



Du {Dl + u 2 ci ) + 9i \ -V, 



d n d 

IT + D *1T 
oy ox 



Viix — D ti [D ti Gn x + oj C i (Guy — Cny)] 



(A4) 



+ { D ti 



ft(l-7) + 4£ 



9 ] „ 



where c si = ( / jPio/'m i n i0 ) 1 ^ 2 is the ion sound velocity. To obtain the second equation 
expressing V • Vji through vn x , we apply the operator D ti d/dy to Equation (A2), D ti d/dz 
to Equation (A3) and add the resulting equations. Then using Equation (24), we find 



! " : 1 <)>r + dz 2 ) 



Dl - d | — + 



V-Vii 



D ti 
+ 



d d 

' / d d\ ( d 2 d 2 \ 

D ti D ti - UciTT ) ~9i [ + ) 

\ ox ay J \ oy 2 oz l J 



(A5) 



Vilx- 



We introduce the following notations: 
K x = 



D tt (D 2 u + u 2 l )+g l (u Jci - + D u — 



D, 



* (1-7) + 4 A 



+ u ci c si — 



2 d_ 
dy j ' 



Mn x — Du [DtiGii x + oj C i (Guy — Cily)] . 



Then Equation (A4) takes the form 



(A6) 



K x v iXx =M iXx + K 2 V ■ v a . 



(AT) 



We further apply the operator K\ to Equation (A5) and use Equation (A7). As a result, 
we find equation for V • Vji 



= K x D tl 
+ 



d d 



+ 



r _ . d <) \ ( if if 

Dti \ D ti — UJ ci — —Qi I -—r + — —z 

ox oy J \oy z oz 2 



. . . d d\ ( d 2 d 2 

Dti[ Du dx~- Uc %) ~ 9i {df + dz^ 



u 81 1 dy 2 dz 2 



M llx (A8) 
\ V-Vii. 
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If we consider Equation (A8) without electromagnetic forces and the background 
magnetic field, we will obtain equation 

HiV ■ v a = 0, 



where 
Hi 



~ D ti + D ti C si 



dy- 



,'dr 2 

+D 2 / 



dx 19i ) dx 



+ 



d 



d 2 



+ 



d 2 



dz 2 dx 



7 - 1) 9i + 



dc 2 - 
dx 



dy 2 



dz 2 



This equation describe the ion sound and internal (ion) gravity waves. When obtaining 
Equation (A8), we have excluded vn x . We also could exclude V • Va- In this case, the term 
dc 2 Jdx in the last term in the expression for ifj would have the sign — . It is connected 
with the form of dependence of va x from V • Vji and vice versa in Equations (A4) and (A5). 



Let us introduce the following notations: 



dz ay 



N ilx = KtDti 
K 3 = 

Then Equation (A8) is given by 



+ 



d 
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Dti Dfi— <-rf r i—— 



n , n 9 9 

dx ay 



~9i 



d^_ 

dy 2 dz 2 



dx 
K 2 -K, 



dy 



-9i 



D 2 - c 2 



dy 2 dz 2 

d 2 d 2 
+ 



dy 2 dz 1 



K 3 V ■ v a = - N ilx . 



(A10) 



Using notations (A6), we can represent the operator K 3 and the value Nn x defined by 
Equations (A9) in the form 



K 3 = uj 2 ci D ti (c 



+Dl 



81 dz 2 

d 2 



-D 2 tl ) +u cl D 2 tt 
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dy 
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d 2 
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dx 



(All) 
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and 



-N, 



ilx 



+ 
+ 



D u D ti 



d_ 

dx 



COr 



d 



d_ 

'dy 
d \ 



d 2 d 2 



dy 2 dz 2 

a 2 



D '-{ D '% + ^'1£) +S '{ D « 
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(A12) 



(Gily — Cily) 



We further assume the following simplifications ~ dc 2 si /dx): 



, , d 2 d 2 d 2 \ - d 



' \ :).,■'• ;). > ;)..2 / ' **' 



dy 2 dz 2 dx 
n 9i d 2 ;,r d 2 



(A13) 



Now, D ti — d/dt + v i0x d/dx. Then the operator A" 3 takes the simple form (see Eq. [All]) 



K 3 = u^Dti ( c 



Sl dz 2 



D u +u a D^ 
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(A14) 



We take into account the last small term on the right hand-side of Equation (A14) to obtain 



some additional terms proportional to co ci 2 in expressions for components of v,i 



A.l. The velocity vn x 



From Equations (A4), (A10), (A12), and (A14), we find the velocity vn X) using 
conditions (A13), 



Vilx = 



( „2 d 2 i „2 d 2 n 2\ 
D ti { C si 9y 2+ C si 9z 2 U ti) I 

-tJ-jlzH \y*i\y ~ ^ily) 



( 



D ti 



-.2 d 2 

1 



D ti) 



(A15) 
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When obtaining solution (A15), we have used some additional conditions 

4|^ + 4|^-a0 » —4^, (Ai6) 



' dy 2 Sl dz 2 1 ) oo C i m dxdy 



( _^U » 9! d 2 D tl (d 2 d 2 d 2 \Dj 
Ul 81 dx) dy ^ D ti u ci dz 2 ' co a Sl \dx 2 ^ dy 2 ^ dz 2 J ' u cl ' 

9i ( d \ ~ 1 f 9 x J 



ujd \oy J \az 

We note that in general expressions, we do not use the local approximation. 



A. 2. The velocity v., 



The velocity vn y , we find from Equation (Al). Applying the operator Du to this 
equation and using Equations (23) and (24), we obtain 

d 



-DtiGu x + D ti u! C iVu v + 



dx 



ft (1"7) + 4^ 



V-v a . 



(A17) 



From Equations (A10), (A12), (A14), (A15), and (A17), we find, using conditions (A13) 
and (A16), 
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Du 



ft (1-7) + 4^ 



(7-2)ft + ^' 



5 2 



G 
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A. 3. The velocity t> il2 



The velocity vn z , we find from Equation (A3). Applying the operator D ti and using 
Equation (24), we obtain 

d d 
D 2 ti v ilz = - 9ifcVnx + DtfQiu + c«i^V • v a . (A19) 

From Equations (A10), (A12), (A14), (A15), and (A19), we find, using conditions given 
above, 
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Note that the small term proportional to ioJ~Gn z in Equation (A20) has appeared due to 
the small term in Equation (A14). 
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B. SOLUTION OF EQUATION (27) 



The components of Equation (27) are the following: 



1 dp el n el 

0= rriig e h G elx + m e uj ce v el y, (Bl ) 

n e0 Ox n e0 

0= + G e iy — C e \y — m e oj ce v e i x , (B2) 

0=-—^ + G elz . (B3) 
n e0 oz 



Below, we find the components of v e ±. 



B.l. The velocity v e \ x 

The velocity v e \ x can be easily found from Equations (B2) and (B3). Differentiating 
Equation (B2) over d/dz and (B3) over d/dy and subtracting one equation from another, 
we obtain 

dv elx d . , d 



m P uj, 



e^ce 



{Gely — G e iy) — -£—G e i z . (B4) 



oz oz ay 



B.2. The value V • v ei 



We can find the value V • v e i from Equation (B2). Applying to this equation operator 
D te and using Equations (36) and (B4), we obtain 
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B.3. The velocity v e i y 



We calculate the velocity v e \ y from Equation (Bl), applying the operator D te and using 
Equations (28) and (36). We further insert in the equation obtained the value V • v el . The 
important point is to differentiate carefully the background electron number density and 
pressure. Then we use Equation (B4) for v e \ x . As a result of calculations, we obtain 



mJjJrj-V 



e UJ ce U e ly 



— Gelx + 



elx i ■ 
C 

^se 



dc 2 

(7-l)<?e+ - 



rrii 



dx \dz J 62 cL 



dx 

(7-1)0 



{jD te + Q) m e u, 
lD te 



(Gely ~ C e iy) (B6) 



+ 



dc 2 



9e 



+mi9e 



dc 2 

dT e0 i n 



( 7 Ae + dx 
to,- d 



(7Ae + m e u ce dy \dz 



7T~ ) G el2 



TeodxBo^Dte + Q) \dz 



l\ B 



It- 



B.4. The velocity v e i 2 



We find this velocity from V • v e i 

dv elz 



= V-v 



dv elx dv, 



dz ' el dx dy 
Using solutions (B4)-(B6), we obtain from Equation (B7) 
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(B8) 



c 2 el midx B Q (D te + Q) \dz / 
The velocity c se i is defined by expression (49). When obtaining Equation (B8), we have 



used condition D te ^> g e k y /uj C i. In this case, D te = d/dt + v e o x d/dx. 



